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Abstract. We consider two mechanisms of spin relaxation in disordered graphene. 
i) Spin relaxation due to curvature spin orbit coupling caused by ripples, ii) Spin 
relaxation due to the interaction of the electronic spin with localized magnetic 
moments at the edges. We obtain analytical expressions for the spin relaxation 
times Tso and tj due to both mechanisms and estimate their values for realis- 
tic parameters of graphene samples. We obtain that spin relaxation originating 
from these mechanisms is very weak and spin coherence is expected in disordered 
graphene up to samples of length L ~ 

1 Introduction 

Two dimensional (2D) graphene, a stable atomic layer of carbon atoms, remained for long 
ellusive among the known crystalline structures of carbon. Only recently, the experimental 
realization of stable, highly crystalline, single layer samples of graphene have been possible 
[112] . Such experimental developments have generated a renewed interest in the field of two 
dimensional mesoscopic systems. The peculiar electronic properties of the hexagonal lattice of 
carbon atoms makes graphene quite different from standard 2D semiconducting heterostructures 
samples. 

The unit cell of graphene is described by two inequivalent triangular sublattices A and 
B intercalated and there are two independent fc-points, K and if', corresponding to the two 
inequivalent corners of the Brillouin zone in reciprocal space. The Fermi level is located at these 
K and K' points and crosses the tt bands of graphene. These two features provide an exotic 
fourfold degeneracy of the low energy (spin-degenerate) states of graphene. These states can be 
described by two sets of two-dimensional chiral spinors which follow the massless Dirac-Weyl 
equation in terms of a pseudo-spin degree of freedom which refers to the fourfold degeneracy 
mentioned. This Dirac-Weyl equation describes the electronic states of the system near the K 
and K' points where the Fermi level is located. At present, there is a large activity in the study 
of the dynamics of this pseudo-spin degree of freedom 

Less attention has been given to the spin so far. The main interactions that could affect 
the spin degree of freedom in graphene seem to be the spin-orbit coupling and exchange in- 
teraction. It is not known to which extent magnetic impurities are present in actual graphene 
samples. Their effect seem small though, as noticed recently when investigating weak localiza- 
tion and universal conductance fluctuations in graphene^. Spin-orbit interaction in graphene 
is supposed to be weak, due to the low atomic number Z = 6 oi carbon. We have investigated 
recently the physics of the spin-orbit interaction in graphene in detail fS]. The effect of other 
interactions as the exchange interaction due to localized moments at the edges of the sample, 
has been discussed elsewhere |6l7l8] . 
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Table 1. Dependence on band structure parameters, curvature, and electric field of the spin orbit 
couplings discussed in the text in the limit Vi <^ V2 (widely separated a bands). The parameters used 

are A ^ 0.264A, £ « 50V/300nm, A = 12meV,"l4pa ~ 4.2eV, 3.63eV, Vpp^ ~ 5.38eV and 

Vpp^ ~ -2.24eV, Vi = 2.47eV, V2 = 6.33eV, a = 1.42A, « ~ 100 A, /i ~ 10 A and i? ~ 50 - lOOnm. 
See Ref. [S] and references therein for details. 



Intrinsic coupling: Ant 1^ {^Y O.OIK 

Rashba coupling (electric field £ « 50y/300nm): As 0.07K 
Curvature coupling: A,ur. ^(v^^„^-v^,^) ^_a_ ^ ^v^y 



2 Spin Orbit 

For the spin-orbit coupling in graphene we have obtained an effective Hamiltonian for the 
TT bands, by second order perturbation theory, which is formally the same as the effective 
Hamiltonian obtained previously from group theoretical methods by Kane and Mele[9]. We find 
for the coupling parameters charaterizing such a Hamiltonian that: i) The intrinsic interaction 
^int ^ lOmK is two orders of magnitude smaller than what was recently estimated by Kane 
and Mele (9i5j . Moreover, we find that for typical values of the electric field as e.g. used by Kane 
and Mele[9], the Rashba- type spin-orbit interaction due to an applied perpendicular electric 
field is A£ ^ 70mK. So spin-orbit coupling for fiat graphene is rather weak. Graphene samples 
seem to have an undulating surface [4]. We have obtained that local curvature of the graphene 
plane induces another type of spin orbit coupling Zicurv Our estimate for the typical observed 
ripples indicates that Z\curv ~ 0.2K. It seems that curvature effects on the scale of the distance 
between neighbouring atoms could increase the strength of the spin-orbit coupling at least one 
order of magnitude with respect to that obtained for a fiat surface. More importantly, this type 
of "intrinsic" coupling will be present in graphene as long as its surface is corrugated even if 
S =0 when As =0. 



2.1 Numerical estimate for Curvature Spin-Orbit. 



"^curv due to the corrugated nature of the samples, is the most important of the spin-orbit 
couplings. The ripples observed seem to be of several A height and a few tens nm laterally. It 
seems possible to identify ripples of lateral size ranging 50nm -lOOnm in Ref. [4 . From these 
ripples we estimate Z\curv 0.2 K[5]. We have checked that this estimate, is in good agreement 
with recent data obtained for the band structure of the tt- and cr-bands [lOj . from which we 
obtain a value Zicurv ~ 0.15K. 

The effective Hamiltonian for graphene, including the spin-orbit coupling reads [5]: 



Tir = y d^rtf-t ^-invpla-^dx + T^aydy] + Aint[f z^zh] + ^[^xSy + fzd-yS^]^ if-, (1) 

where hvp — \/37oa/2, a ^ 2.46 A being the lattice constant for graphene, 70 ~ 3eV the 
McClure intralayer coupling constant, Z\int ~ O.OIK, Aj^ ~ As + Z\curv the Rashba-Curvature 
couphng (RCC), where As ^ 0.07K for £ « 50V/300nm and Z\curv ~ 0.2K. Tabled summarizes 
the main results obtained for the effective spin-orbit couplings in Ref. [5]. 
In the following, we consider £ — 0, An = Acuiv and Z\int ^ ^i?: 

Hr = J d^.^.) {-^nMoA + f.^ydy] + A,[,.sy ^ ^zOyS.]) ^(x) (2) 
We re- write TLt = Ho + Ti^ where: 

Ho = y" d^x^\^) (^-ihvFi^Jx + Tz^ydyij (x), (3) 



Will be inserted by the editor 



3 



and 

HRix)^ J d'xWH^){ARi^)[a,Sy+T,ayS,])Wix.) = J d^xAR{x.)^\^)rW{^). (4) 

The coupling Ar('x.) depends on the position of each ripple, as the corrugations depend on the 
position and the defined curvature may change sign across the sample randomly. The average 
radius of curvature R is zero and Aji{'x.) can be considered as a quenched Gaussian variable: 

(Z\fl(x)) = 0; (z\fl(x)Z\fl(x )) = AlS\jc-x')f, (5) 

where I provides the scale for both the average size of each ripple and the average distance 
between ripples and < > denotes ensemble averaging over different curved samples. 

Following [6j the self-energy due to this type of "topological disorder" reads in the SCBA: 

Sr = (AR{q)rG{p + q, e)rAR{-q)) = J ^Zi^(q)f G(p + q, 6)f Zifl(-q). (6) 

Eq.(I5]) implies in momentum space: 

{AR{<DARi-<i))^Alf, (7) 

and finally we have: 

Er = AlPTr[ J ^^rG(k,e)r], (8) 

where Tr[ ] denotes trace over all possible spin and pseudo-spin indexes and G(k, e) is the 
Green's function averaged over disorder. 

We are interested in the components of the self-energy matrix which describe spin flip 
processes: 

Im.[s+Srs-] - All^h+Im[ J ^AA_Tr[a+G(k, (9) 

and 

Im[s.SrH] - All^h^Iml J ^^Tr[a^G{k, e)a^]]S+, (10) 

where now Tr[ ] is only trace over pseudo-spin indexes. Note that the processes that provide 
spin flip, Eq.® and Eq.([TO]), are contributions in which K and K' are mixed by disorder. 
The only finite contributions to spin flip are contributions e.g. where = +1 for the flrst 
r and Tz = —1 for the second F in Eq.®. A process t^i requires K ^ K' and for 
requires K' K. Summing up both contributions, we finally get for the imaginary part of 
Scurv = s+SpS- + s_Z'r.s+: 

/m[i:,„™] = Alphm[ I J^Tr[G{k, e)] ] = AUMe)- (H) 

Note that there is no extra valley degeneracy and in this case spin and valley degeneracy are 
coupled. This means that i^(e) in Eq. (fTT|) is 1/2 of the total density of states for a disordered 
graphene sample ^'r(e) = 5s5t)|e|/(27r(?iWF)^), where gs = gv = are the spin and valley 
degeneracies respectively. At zero applied perpendicular electric field, Ar = Zicurv and we 
obtain the spin relaxation time due to curvature spin-orbit coupling as: 

— = Im[s+Srs-] = Al^„fv{e). (12) 

TSO 

In clean graphene the "spin-valley" density of states for energies close to the Dirac point, e ~ 0, 
reads: 
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The presence of vacancies or dislocations leads to a finite density states at the Dirac point, 
e = [8], which in our "spin- valley" case reads: 

.(. = 0).^. (14) 

1 is the mean free path associated with the dislocations or vancancies in the system |8|. We 
assume that due to dislocations or vacancies, there is finite intervalley scattering in the sam- 
ple at length scales of the order of 1. Also if the graphene sample is tightly coupled to the 
insulating Si02 substrate, this generates atomically sharp scatterers which lead to intervalley 
scattering [IjJ. Using hvp = fia ^ b.ieVA, where t is the hopping between neighboring atoms 
in the graphene lattice, 1 = lO^a, where a is the lattice constant of graphene, a = 2.46A, 
I = lOOA and Z\curv ~ Q.2K ^ 1.742 x lQ~^eV , we finally have for the spin-relaxation time 
close to the Dirac points K : 

^so = - - 0.282 X 10-S. (15) 

curv^ 



3 Localized moments 

Lattice distorsions and sample edges in graphene can induce localized states at the Fermi 
energy, leading to the existence of local moments [7^. The RKKY interaction between these mo- 
ments is always ferromagnetic due to the semimetallic properties of graphene. Instability of the 
paramagnetic phase due to electron-electron interaction has been also proposed as mechanisms 
leading to ferromagnetism in graphenejiSj. Also, spin-flip for spin splitted edge channels in the 
QH regime has been estimated to be rather weak|12|. 

The interaction between the electronic spin and the localized magnetic moments is given by 
Heisenberg Hamiltonian: 

Tij^ JS-s (16) 

where 

J=^^C/ (17) 

where t ^ U because for graphene /fivp = {e'^)/{Wa) = U/W ~ 1. 

It is straightforward to obtain the spin-relaxation time due to the interaction between the 
electronic spins and the localized moments: 

A. 5(5 + 1),^^. (18) 

N denotes the total probability for an electronic quasiparticle to reach the edges and interact 
with the localized moment. Using N = N fia/L)^ , where N is the fraction of localized moments 
per edge and / is the fraction of edge which is a zigzag where localized states occur, we finally 
have: 

"•^ = K0 ir,s{sli)Nf (^) ■ ^^^^ 

Taking N — 0.3, i.e. there is a localized moment every three atoms of a zizag edge, / = 0.2 as 
the average fraction of a sample edge which is of zigzag type and for the sample size L = 100/im, 
we have 

Tj - 0.16 X lO^^s. (20) 
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4 DifFusion coefficient in graphene 

From the mobility of the graphene samples /i — a /(en), n being the electron(hole) density and e 
the electronic charge, and using the Einstein relation a — e^h'{ep)D, we obtain for the diffusion 
coefficient: 

D^t^^. (2f) 

e 

For /i ^ 10'*cTO^/(Fs)[T], n = lO^^cm"^ [13], we obtain a very rough estimate for the diffusion 
coefficient: 

D-2xfO"Ws"\ (22) 

Under the assumption that transport in disordered graphene is diffusive, we obtain the following 
order of magnitude for the spin-flip lengths: 

Aso = VDtso Aj = ^/D^J - VO.2 IQ-^ x 2 IQ-^ - 2fim. (23) 



5 Conclusions 

We have investigated two mechanisms of spin relaxation in disordered graphene: i) Spin-flip 
due to the spin-orbit interaction caused by ripples in the surface of graphene and ii) Spin-flip 
due to the interaction of the spin with magnetic moments localized at the sample edges. We 
have found out that spin-relaxation from these mechanisms is very weak with a relaxation time 
TsoiJ) ^ 0.2 X 10~^s. By estimating the diffusion coefficient in graphene samples, we obtain 
that the spin-flip length due to spin-orbit and/or interaction with localized moments Xso(J) is 
of the orther of ~ f/im. So spin relaxation should start being important if the sample size L 
becomes of the other of L '--^ l^m. This shows that spin coherence in graphene can be preserved 
for long distances. These results support perspectives for spintronics applications of graphene. 
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grants no 162742/vOO, 1585f8f/43f and f 158547/43f . F. G. acknowledges funding from MEC 
(Spain) through grant FIS2005-05478-C02-01 and the European Union Contract 12881 (NEST) 
and the Comunidad de Madrid, through the program CITECNOMIK, CM2006-S-0505-ESP- 
0337. 
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